Surface location of sodium atoms attached to 3 He nanodroplets 
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We have experimentally studied the electronic 3p «— 3s excitation of Na atoms attached to 3 He 
droplets by means of laser-induced fluorescence as well as beam depletion spectroscopy. From the 
similarities of the spectra (width/shift of absorption lines) with these of Na on 4 He droplets, we 
conclude that sodium atoms reside in a "dimple" on the droplet surface and that superfluid-related 
effects are negligible. The experimental results are supported by Density Functional calculations 
at zero temperature, which confirm the surface location of sodium on 3 He droplets, and provide a 
detailed description of the "dimple" structure. The calculated shift of the excitation spectra for the 
two isotopes is in good agreement with the experimental data. 
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Detection of laser-induced fluorescence (LIF) and 
beam depletion (BD) signals upon laser excitation pro- 
vides a sensitive spectroscopic technique to investigate 
electronic transitions of chromophores attached to 4 He 
nanodroplets Q. While most of atomic and molecu- 
lar dopants migrate to the center of the droplet, alkali 
atoms (and alkaline earth atoms to some extent Q) have 
been found to reside on the surface of 4 He droplets, as 
evidenced by the much narrower and less shifted spec- 
tra when compared to those found in bulk liquid 4 He 
LH IE 0- This result has been confirmed by Den- 
sity Functional (DF) Q and Path fntegral Monte Carlo 
(PIMC) |8j] calculations, which predict surface binding 
energies of a few Kelvin, in agreement with the measure- 
ments of detachment energy thresholds using the free 
atomic emissions |9|- The surface of liquid He is only 
slightly perturbed by the presence of the impurity, which 
produces a "dimple" on the underlying liquid. The study 
of these states can thus provide useful information on 
surface properties of He nanodroplets complementary to 
that supplied by molecular-beam scattering experiments 

Although the largest amount of work has been devoted 
to the study of pure and doped 4 He nanodroplets -see 
and Refs. therein-, the only neutral Fermi systems 
capable of being observed as bulk liquid and droplets are 
made of 3 He atoms, and for this reason they have also at- 
tracted the interest of experimentalists and theoreticians 
ElQEEHEIllil. We recall that while 4 He 
droplets, which are detected at an experimental temper- 
ature (T) of ~ 0.38 K, are superfluid, these containing 
only 3 He atoms, even though detected at a lower T of ~ 
0.15 K, do not exhibit superfluidity pi) . 

The behavior of molecules in He clusters is especially 
appealing. In particular, probes at the surface of the 
droplets are desirable because they allow to investigate 
the liquid-vacuum interface as well as droplet surface 
excitations. The latter are of interest in the compari- 



son of the superfluid vs. normal fluid behavior, particu- 
larly because the Bose-Einstein condensate fraction has 
been calculated to approach 100% on the surface of 4 He 
[22| . Small 3 He drops are difficult to detect since, as a 
consequence of the large zero-point motion, a minimum 
number of atoms is needed to produce a selfbound drop 
0130. Microscopic calculations of 3 He droplets are 
scarce, and only concern the ground state (GS) struc- 
ture (lH l2(il ] . GS properties and collective excitations of 
3 Hc droplets doped with some inert atoms and molecular 
impurities have been addressed within the Finite Range 
Density Functional (FRDF) theory that has proven 
to be a valuable alternative to Monte Carlo methods, 
which are notoriously difficult to apply to Fermi systems. 
Indeed, a quite accurate description of the properties of 
inhomogeneous liquid 4 He at T — has been obtained 
within DF theory p3j , an d a similar approach has fol- 
lowed for 3 He (see [ljj,|2^ and Refs. therein). 

The experiments we report have been performed in 
a helium droplet machine used earlier for LIF and BD 
studies, and is described elsewhere |2|. Briefly, helium 
gas is expanded under supersonic conditions from a cold 
nozzle forming a beam of droplets traveling freely under 
high vacuum conditions. The droplets are doped down- 
stream employing the pick-up technique: in a heated 
scattering cell, bulk sodium is evaporated in such a way 
that, on average, a single metal atom is carried by each 
droplet. LIF absorption spectra of doped droplets are 
recorded upon electronic excitation using a continuous 
wave ring-dye laser and detection in a photo multiplier 
tube (PMT). Since electronic excitation of alkali-doped 
helium droplets is eventually followed by desorption of 
the chromophorc, BD spectra can be registered by a 
Langmuir- Taylor surface ionization detector [25| • Phase- 
sensitive detection with respect to the chopped laser or 
droplet beam was used. For that reason the BD sig- 
nal (cf. Fig. |2J, i.e. a decrease in intensity, is directly 
recorded as a positive yield. For these experiments, a 
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FIG. 1: Laser-induced fluorescence signal as a function of noz- 
zle temperature forming 3 He droplets (black) in comparison 
to 4 He droplets (grey). In both runs a stagnation pressure of 
20 bar was used; nozzle diameter was 5 /im. Normalization is 
such that the plot gives the correct relative intensities. 
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FIG. 2: Beam depletion spectra of Na atoms attached to 
3 He/ 4 He nanodroplets. The vertical lines indicate the po- 
sitions of the two components of the Na gas-phase 3p <— 3s 
transition. 



new droplet source was built to provide the necessary 
lower nozzle temperatures to condense 3 He droplets. Ex- 
panding Pq = 20 bar of helium gas through a nozzle 5 /zm 
in diameter, we now can establish temperatures down to 
7.5 K using a two-stage closed cycle refrigerator (Sumit- 
omo Heavy Industries, Model: RDK-408D). In this way, 
without needing any liquid helium or nitrogen for pre- 
cooling or cold shields, stable beam conditions can be 
utilized over several days. 

Fig. compares the number of Na-doped 3 He vs. 4 He 
droplets expanding 20 bar helium as a function of the noz- 
zle temperature To. The absolute number densities in the 
maxima of both distributions are quite similar. The for- 
mation of droplets from a supersonic expansion is well 
known and thoroughly discussed in the literature: Be- 
sides the determination of absolute cluster sizes 0, |2(| , 
the size dependence as a function of source conditions 
has already been studied 0, ■ The low temperature 
cutoff appears at source conditions where the isentrope 
of the expansion hits the helium critical point [2^; the 
disappearance at high temperatures just means that the 
droplets are getting to small to carry the dopant. Here, 
the critical size is determined by the thermal energy dur- 
ing pick-up which leads to evaporation of helium atoms 
and a destruction of small droplets. For the spectro- 
scopic measurements presented in the following, we set 
T = 11 K for 3 He, and T = 15 K for 4 He. These condi- 
tions are expected to result in comparable mean cluster 
sizes around 5000 atoms per droplet 0, |2(j . As far as 
our results are concerned, Fig. ^ demonstrates that we 



see the correct formation as well as droplet sizes of the 
different isotopes and that we use suitable source condi- 
tions to guarantee comparable sizes. 

In Fig.|21the absorption spectrum of Na atoms attached 
to 3 He nanodroplets is shown in comparison to Na-doped 
4 He droplets. We present here the BD spectra because 
they do not contain the strong fluorescence background 
lines of free Na atoms which cover the crucial steep in- 
crease of the droplet spectrum. Moreover, LIF does not 
always represent the total absorption spectrum because 
it relies on the emission of a photon in the spectral range 
of the PMT. Hence, absorption processes followed by ei- 
ther radiationless decay or emission of photons in the 
infrared spectral region are suppressed. The latter has 
been observed in LIF spectra where alkali-helium exci- 
plexes form upon excitation of alkali atoms on the surface 
of 4 He droplets UEH 

In our experiment Na 3 He exci- 
plexes are formed in the same way as their Na 4 He coun- 
terparts. This became immediately obvious because we 
were able to discriminate the corresponding red-shifted 
emission intensities. However, as far as the measured 
absorption spectra of Na@ 3 Hejv are concerned, the LIF 
data are very well in accord with the BD absorption. 

The outcome of the spectrum of Na attached to 3 He 
nanodroplets is very similar to the spectrum on 4 Hc 
droplets. The asymmetrically broadened line is almost 
unshifted with respect to the gas-phase absorption. This 
absence of a shift immediately confirms the surface lo- 
cation because atoms embedded in bulk superfluid he- 
lium are known to evolve large blue-shifts of the order 
of a couple of hundreds of wavenumbers and much more 




FIG. 3: Equidensity lines in the x — z plane showing the stable 
state of a Na atom (cross) on a He2ooo droplet. The 9 inner 
lines correspond to densities 0.9po to O.lpo, and the 3 outer 
lines to 10~ 2 p , 10~ 3 po, and 10~ 4 p {pa = 0.0163 A' 3 for 
3 He, and 0.0218 A" 3 for 4 He). 

broadened absorption lines ,30]. A blue shift is a conse- 
quence of the repulsion of the helium environment against 
the spatially enlarged electronic distribution of the ex- 
cited state ("bubble effect"). The interaction towards the 
3 He droplets appears to be slightly enhanced, evidenced 
by the small extra blue shift of the spectrum compared 
to the 4 He spectrum. In a simple picture this means 
that more helium atoms are contributing or, in other 
words, a more prominent "dimple" interacts with the 
chromophore. The upper halves of the spectra are almost 
identical, when shifting the 3 He spectrum by 7.5 ± 1 cm -1 
to lower frequencies. The influence of different droplet 
sizes does not affect this observation: A dependence of 
the 4 He spectra varying the droplet size is already shown 
in Q. Probing smaller droplets narrows the spectrum 
but does not shift the observed doublet structure. The 
shift of the 3 He spectrum would be even more conspic- 
uous when comparing smaller droplets. Unfortunately, 
the corresponding small 3 He droplet sizes are experimen- 
tally not accessible because within the plateau region of 
the droplet yield shown in Fig. ^ the droplet sizes are 
almost unchanged [ll|. The increase in width (FWHM) 
of the 3 He spectrum is only 7%. Taking into account the 
just mentioned droplet size dependence of the width, this 
difference might not even be significant. Regarding the 
substructure of the line, the only notable exception is the 
absence of the red-shifted shoulder, which is observed in 
the case of 4 He and marked with an arrow in Fig. [3 This 
feature, which is even more pronounced in the absorption 
of Li-doped 4 He droplets Q has not been interpreted yet. 
The shift with respect to the maximum of the absorption 
line is « 20 cm -1 , too high in energy to be attributed to 
an excited compressional or surface mode of the droplet 
at 0.38 K [H|33. 

FRDF calculations at T = confirm the picture emerg- 
ing from the measurements, i.e. the surface location of Na 
on 3 He nanodroplets causing a more pronounced "dim- 
ple" than in 4 He droplets. We have investigated the sta- 
ble configurations of a sodium atom on both 3 He and 
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FIG. 4: Density profiles along a line connecting the impurity 
to the center of the cluster showing the equilibrium configu- 
ration of a Na atom (cross) on He2ooo nanodroplets. 

4 He clusters of different sizes. The FRDF's used for 3 He 
and 4 He are described in [32I l33l | . The large number of 
3 He atoms we are considering allows to use the extended 
Thomas- Fermi approximation |16| . The minimization of 
the energy DF's with respect to density variations, sub- 
ject to the constraint of a given number of He atoms N, 
leads to Euler-Lagrange equations whose solution give 
the equilibrium particle densities p(v). These equations 
have been solved as indicated in The presence of 

the foreign impurity is modeled by a suitable potential 
obtained by folding the helium density with a Na-He pair 
potential. We have used the potential proposed by Patil 
[35| to describe the impurity- He interactions. Potential 
energy curves describing the He-alkali interaction have 
been calculated recently by ab-initio methods |||, and 
found to agree very well with the Patil potential. 

Fig. 13 shows the equilibrium configuration for a Na 
atom adsorbed onto He2ooo clusters. For a given TV, the 
size of the 3 He at droplet is larger than that of the 4 Hejv 
droplet -an obvious consequence of the smaller 3 He satu- 
ration density-. Comparison with the stable state on the 
4 He2ooo cluster shows that, in agreement with the experi- 
mental findings presented before, the "dimple" structure 
is more pronounced in the case of 3 He, and that the Na 
impurity lies inside the surface region for 3 He and out- 
side the surface region for 4 He (we recall that the surface 
region is usually defined as that comprised between the 
radii at which p = O.lpo and p = 0.9po, where po is the 
He saturation density jloL fill flq|). We attribute this to 
the lower surface tension of 3 He (0.113 K/A 2 ) as com- 
pared to that of 4 He (0.274 K/A 2 ), which also makes the 
surface thickness of bulk liquid and droplets larger for 
3 He than for 4 He 0, ^] . The Na-droplet equilibrium 
distance, here defined as the 'radial' distance between the 
impurity and the point where the density of the pure drop 
would be p ~ po/2, is R ~ 1.1 A for 3 He, and R ~ 3.6 
A for 4 He. For the larger droplets we have studied, R is 
nearly N- independent. A related quantity is the defor- 
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mation of the surface upon Na adsorption, which can be 
characterized @ by the "dimple" depth, £, defined as the 
difference between the position of the dividing surface at 
P ~ Po/2, with and without impurity, respectively. We 
find £ ~ 4.5 Afor 3 He, and f ~ 2.1 A for 4 He. Fig. H 
shows the density profiles for Na@He2ooo- Note the more 
diffuse liquid-vacuum interface for the 3 He droplet far 
from the impurity, and the occurrence of more marked 
density oscillations (with respect to the 4 He case) where 
the softer 3 He surface is compressed by the adsorbed Na 
atom. Our calculations thus yield the detailed structure 
of the liquid around the impurity, which is an essential 
ingredient for any line-shift calculation of the main elec- 
tronic transitions in the adatom Q and for the under- 
standing of dynamical processes which already have beer 
observed in time-dependent experiments 36]. 

We have obtained the shift between the 3 He and 4 Hc 
spectra in Fig. [21 within the Frank-Condon approxima- 
tion, i.e. assuming that the "dimple" shape does nol 
change during the Na excitation The shift is calculatec 
within the model given in |37|. evaluating Eqs. 5 anc 
6 therein, both for 3 He and 4 He. We used the excitec 
state A 2 n and B 2 S potentials of |8j because their Na- 
He GS potential is very similar to the Patil potential we 
have used to obtain the equilibrium configurations. Foi 
our largest droplet (N = 2000) we find that the 3 He 
spectrum is blue-shifted with respect to the 4 He one bj 
6.4 cm -1 , in good agreement with the experimental value 
of 7.5 ± 1 cm -1 as extracted from Fig. [2 

Solvation energies Sno, — -E(Na@HeAr) — E^Hejv) an 
shown in Fig. [5] as a function of droplet sizes. More neg- 
ative values of Sno, indicate increased binding of the Ne 
atom to the droplet surface. The lines are a fit of the 
form 

S Na (N) = So + ^3 + ^ (1) 

with So = -12.1(-12.5), Si = -28.3(-31.6), and S 2 = 
37.3(37.6) K for 4 He( 3 He). To compare with "exact" 
PIMC result 0, we have calculated Na@ 4 He 300 . Part 
of the small difference between the two values has to be 
attributed to our neglecting of the Na zero-point energy. 
The value Sjv a ~ —12 K has been obtained within FRDF 



theory for Na adsorbed on the planar surface of 4 He , 
which corresponds to the N = oo limit in this figure. 

Our results show that Na adsorption on 3 He droplets 
occurs in very much the same way as in the case of 
4 He, i.e., the adatom is located on the surface, though 
in a slightly more pronounced "dimple" . The similari- 
ties in the experimental spectra are certainly remarkable 
for two apparently very different fluids, one normal and 
the other superfluid, and clearly indicate that superflu- 
idity does not play any substantial role in the processes 
described here. This is likely a consequence of the very 
fast time scale characterizing the Na electronic excita- 
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FIG. 5: Na solvation energy as a function of N. The cross is 
the PIMC result of 0. 

is due to the different structure of the "dimple", which 
accounts for the small shift in their spectra observed in 
the experiments and found in our calculations as well. 

We thank Flavio Toigo for useful comments. This 
work has been supported by grants MIUR-COFIN 
2001 (Italy), BFM2002-01868 from DGI (Spain), and 
2001SGR-00064 from Generalitat of Catalunya as well 
as the DFG (Germany). 



[1] F. Stienkemeier and A. F. Vilesov, J. Chem. Phys. 115, 
10119 (2001). 

[2] F. Stienkemeier, F. Meier, and H. O. Lutz, J. Chem. 

Phys. 107, 10816 (1997); Eur. Phys. J. D 9, 313 (1999). 
[3] F. Stienkemeier, J. Higgins, C. Callegari, S. I. Kanorsky, 

W. E. Ernst, and G. Scoles, Z. Phys. D: At., Mol. Clusters 

38, 253 (1996). 
[4] F. Stienkemeier, W. E. Ernst, J. Higgins, and G. Scoles, 

J. Chem. Phys. 102, 615 (1995); 
[5] C. Callegari J. Higgins, F. Stienkemeier, and G. Scoles, 

J. Phys. Chem. A 102, 95 (1998). 
[6] F. Briihl, R. Trasca, and W. E. Ernst, J. Chem. Phys. 



115, 10220 (2001). 
[7] F. Ancilotto, E. Cheng, M. W. Cole, and F. Toigo, Z. 

Phys. B: Condens. Matter 98, 323 (1995). 
[8] A. Nakayama and K. Yamashita, J. Chem. Phys. 114, 

780 (2001). 

[9] J. Reho, C. Callegari, J. Higgins, W. E. Ernst, K. K. 
Lehmann, and G. Scoles, Faraday Discuss. 108, 161 
(1997). 

[10] F. Dalfovo, J. Harms, and J. P. Toennies, Phys. Rev. B 

58, 3341 (1998). 
[11] J. Harms, P. J. Toennies, M. Barranco, and M. Pi, Phys. 

Rev. B 63, 184513 (2001). 



■5 



[12] J. P. Toennies and A. F. Vilesov, Annu. Rev. Phys. 
Chem. 49, 1 (1998); K. K. Lehmann and G. Scoles, Sci- 
ence 279, 2065 (1998), J. P. Toennies, A. F. Vilesov, and 
K. B. Whaley, Physics Today 54, 31 (2001). 

[13] C. Callegari, A. Conjusteau, I. Reinhard, K. K. 
Lehmann, G. Scoles, and F. Dalfovo, Phys. Rev. Lett. 83, 
5058 (1999); 84, 1848(E) (2000); Y. Kwon, P. Huang, M. 
V. Patel, D. Blume, and K. B. Whaley, J. Chem. Phys. 
113, 6469 (2000); E. W. Draeger and D. M. Ceperley, 
Phys. Rev. Lett. 90, 65301 (2003); R. Briihl, R. Guardi- 
ola, A. Kalinin, O. Kornilov, J. Navarro, T. Savas, and 
J. P. Toennies, ibid. 92, 185301 (2004). 

[14] J. Harms, M. Hartmann, J. P. Toennies, A. F. Vilesov, 
and B. Sartakov, J. Mol. Spectrosc. 185, 204 (1997). 

[15] V. R. Pandharipande, S. C. Pieper, and R. B. Wiringa, 
Phys. Rev. B 34, 4571 (1986). 

[16] S. Stringari and J. Treiner, J. Chem. Phys. 87, 5021 
(1987). 

[17] M. Barranco, J. Navarro, and A. Poves, Phys. Rev. Lett. 

78, 4729 (1997). 
[18] LI. Serra, J. Navarro, M. Barranco, and Nguyen Van Giai, 

Phys. Rev. Lett. 67, 2311 (1991); S. Weisgerber and P.-G. 

Reinhard, Z. Phys. D: At., Mol. Clusters 23, 275 (1992); 

M. Barranco, D. M. Jezek, E. S. Hernandez, J. Navarro, 

and LI. Serra, ibid. 28, 257 (1993). 
[19] F. Garcias, LI. Serra, M. Casas, and M. Barranco, J. 

Chem. Phys. 108, 9102 (1998); ibid. 115, 10154 (2001). 
[20] R. Guardiola, Phys. Rev. B 62, 3416 (2000). 
[21] S. Grebenev, J. P. Toennies, and A. F. Vilesov, Science 

279, 2083 (1998). 
[22] A. Griffin and S. Stringari, Phys. Rev. Lett. 76, 259 

(1996); E. W. Draeger and D. M. Ceperley, ibid. 89, 



015301 (2002). 

[23] F. Dalfovo, A. Lastri, L. Pricaupenko, S. Stringari, and 

J. Treiner, Phys. Rev. B 52, 1193 (1995). 
[24] E. S. Hernandez and J. Navarro, in Microscopic Ap- 
proaches to Quantum Liquids in Confined Geometries, 

E. Krotscheck and J. Navarro Eds. pag. 261 (World Sci., 

Singapore, 2002). 
[25] F. Stienkemeier, M. Wewer, F. Meier, and H. O. Lutz, 

Rev. Sci. Instr. 71, 3480 (2000). 
[26] J. Harms and J. P. Toennies, unpublished results. 
[27] A. Scheidemann, J. P. Toennies, and J. A. Northby, Phys. 

Rev. Lett. 64, 1889 (1990). 
[28] J. Harms, J. P. Toennies, and E. Knuth, J. Chem. Phys. 

106, 3348 (1997) 
[29] J. Reho, J. Higgins, C. Callegari, K. K. Lehmann, and 

G. Scoles, J. Chem. Phys. 113, 9686 (2000). 
[30] Y. Takahashi, K. Sano, T. Kinoshita, and T. Yabuzaki, 

Phys. Rev. Lett. 71, 1035 (1993). 
[31] M. Casas and S. Stringari, J. Low Temp. Phys. 79, 135 

(1990); S. A. Chin and E. Krotschek, Phys. Rev. B 45, 

852 (1992); ibid. 52, 10405 (1995). M. Barranco and E. 

S. Hernandez, ibid. 49, 12078 (1994). 
[32] M. Barranco, M. Pi, E. S. Hernandez, and J. Navarro, 

Phys. Rev. B 56, 8997 (1997). 
[33] R. Mayol, M. Pi, M. Barranco, and F. Dalfovo, Phys. 

Rev. Lett. 87, 145301 (2001). 
[34] M. Barranco, M. Guilleumas, E. S. Hernandez, R. Mayol, 

M. Pi, and L. Szybisz, Phys. Rev. B 68, 024515 (2003). 
[35] S. H. Patil, J. Chem. Phys. 94, 8089 (1991). 
[36] O. Biinermann et al., unpublished (2004). 
[37] S. I. Kanorsky, M. Arndt, R. Dziewior, A. Weis, and T. 

W. Hansen Phys. Rev. B 50, 6296 (1994). 



